
ANALYSIS OF STRUCTURAL RESPONSE 
WITH DIFFERENT FORMS OF DAMPING 


by R. R. Reed 


George G Marshall Space Flight Center 
Huntsville, Ala. 

• WASHINGTON, 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



□13DbS0 


ANALYSIS OF STRUCTURAL RESPONSE WITH 
DIFFERENT FORMS OF DAMPING 
By R. R. Reed 

George C. Marshall Space Flight Center 
Huntsville, Ala. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 




TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

ANALYTICAL REPRESENTATION OF DAMPING 2 

ENERGY LOSSES AND HYSTERESIS CURVES 6 

TRANSIENT RESPONSE 15 

CONCLUSIONS 22 

APPENDIX 41 

REFERENCES 48 


iii 



LIST OF ILLUSTRATIONS 


Figure Title Page 

1. Nonlinear Spring Force 3 

2. Energy Loss Ratio versus Frequency. . 10 

3. Energy Loss Ratio versus Amplitude n 

4. Single Hysteresis Loops 16 

5. Double Hysteresis Loops 17 

6. -Displacement versus Nondimensional Time - Case I 

x 3 

(x + 0. 04x + x - — = 0. 2) 24 

7. Displacement versus Nondimensional Time - Case I 

x 3 

(x + 0. 2x + x - — = 0. 2) 24 

O 

8. Displacement versus Nondimensional Time - Case I 

x 3 

(x + 0. 04x + x - — = 0, x(t = 0) = 0. 2) 25 

o 

9. Displacement versus Nondimensional Time - Case I 

X 3 

(x + 0. 2x + x - — = 0, x(t = 0) = 0. 2) 25 

o 

10. Force versus Displacement • Case I 

X 3 

(x + 0. 2x + x - — = 0. 2) 26 

11. Force versus Displacement - Case I 

x 3 „„ 

(x + 0. 2x + x - — = 0, x( r = 0) = 0. 2) 26 

12. Displacement versus Nondimensional Time - Case II 

x 3 

(x +0.04|x|(sgn x) + x - — = 0.2) . . 27 


iv 



LIST OF ILLUSTRATIONS (Continued) 

Figure Title Page 

13. Displacement versus Nondimensional Time - Case II 

x 3 

(x + 0. 2 | x| ( sgn x) + x - — = 0. 2) 27 

14. Displacement versus Nondimensional Time - Case II 

X 3 oo 

(x + 0. 04| x| ( sgn x) + x - — = 0, x(t = 0) = 0. 2) 28 

15. Displacement versus Nondimensional Time - Case II 


(x + 0. 2 |x| (sgn x) + x - — = 0, x(r = 0) = 0. 2) 28 

16. Force versus Displacement - Case II 

x 3 

(x + 0. 2 1 x| (sgn x) + x - — = 0.2) 29 


17. Force versus Displacement - Case II 

x 3 

(x + 0. 2 |x| (sgn x) + x - — = 0, x(r = 0) = 0. 2) 29 


18. Displacement versus Nondimensional Time - Case III 

( x + 0. 04k 2 + x - ~ = 0. 2) 30 

19. Displacement versus Nondimensional Time - Case III 

(x + 0. 2x? + x - ~ = 0. 2) 30 

O 

20. Displacement versus Nondimensional Time - Case III 

x 3 

(x + 0.04x 2 + x - f - = 0, x(t = 0) = 0. 2) 31 

21. Displacement versus Nondimensional Time - Case III 

X 3 

(x + 0. 2x 2 + x - — = 0, x(t = 0) = 0. 2) 31 

o 

22. Force versus Displacement - Case Ell 

(x + 0.2x 2 + x- ~ = 0.2) 32 

o 


V 



LIST OF ILLUSTRATIONS (Continued) 


Figure 

Title 

Page 

23. 

Force versus Displacement - Case in 

X 3 

(x + 0. 2x 2 + x - — = 0, x(t = 0) = 0. 2) 

d 

. . . ,32 

24. 

Displacement versus Nondimensional Time - Case IV 

X 3 

(x + 0.04jx|x + x - — = 0. 2) 

. . . 33 

25. 

Displacement versus Nondimensional Time - Case IV 

X 3 

(x + 0. 2 |x|x + x - — = 0.2) 

d 

. . . 33 

26. 

Displacement versus Nondimensional Time - Case IV 

X 3 

(x + 0. 04| x| x + x - — = 0, x(t = 0) = 0. 2) 

d 

. . . 34 

27. 

Displacement versus Nondimensional Time - Case IV 

X 3 

(x + 0. 2|x|x + x - — = 0, x(r = 0) = 0. 2) 

o 

. . . 34 

28. 

Force versus Displacement r Case IV 

X 3 

(x + 0. 2 |x|x + x - — = 0. 2) 

o 

. . . 35 

29. 

Force versus Displacement - Case IV 

X 3 

( x + 0. 2 1 x | x + x - — =0, x( r = 0) = 0. 2) 

d 

. . . 35 

30. 

Displacement versus Nondimensional Time - Case V 

X 3 

(x + 0. 0392(1 + 3^)x + x - — = 0. 2) 

d 

. . 36 

31. 

Displacement versus Nondimensional Time - Case V 

X 3 

(x + 0. 10( 1 + x 2 )x + x - — = 0. 2) 

. . 36 


32. Displacement versus Nondimensional Time - Case V 

X 3 

( x + 0. 039;2( i + x 2 ) x + x - — = 0, x( r = 0) =0.2) 37 

u 


vl 



LIST OF ILLUSTRATIONS (Concluded) 


Figure Title Page 

33. . Displacement versus Nondimensional Time - Case V 

x 3 

(x + 0. 10( 1 + x?) x + x - — = 0, x(t = 0) = 0. 2) 37 

u 

34. Force versus Displacement - Case IV 

x 3 

( x + 0. i ( 1 + x 2 ) x + x - — = 0. 2) 38 

o 

f 

35. Force versus Displacement {x + 0. i( i + x*)& + x -~r** 0, 

x (t = 0) = 0. 2 ' 38 

A. i. Phasor Diagrams 47 

LIST OF TABLES 

Table Title Page 

A 

I . Analytical Representation of Damping v 

II. Energy Losses - Sinusoidal Response *7 

III . Energy Loss Ratios 8 

IV. Transient Energy Losses - Free Response 43 

V. Transient Energy Losses - Step Response 14 

VI . Solution Exponent, Frequency, and Steady State 

Amplitude - Step Input 19 

VII. . Solution Exponent and Frequency - Free Oscillation .20 

VIII. Maximum and Minimum x versus r Values - Step Input 21 

IX. Maximum and Minimum x versus r Values - Free 

Oscillation 21 

X. Summary of Curves 39 

vii 



DEFINITION OF SYMBOLS 


Symbol 


a 

A o 

C()j Cj, C2 j C3 

C lr 

K 

M 

P 


R 

r 

T 

W 


a , b 
n n 


n, m 
t 


x„ 


x 10> x 20 
AW 


Definition 

amplitude of the response 

the n* 11 complex coefficient of the exponential series 

steady state response amplitude 

damping coefficients 

a reference value of Cj 

linear spring stiffness coefficient 

generalized mass 

magnitude of step input force 

energy loss ratio 

referenced energy loss ratio 

the initial time reference for energy measurements 

total energy content 

real and imaginary parts of 

whole numbers 

independent variable time 

generalized mass displacement 

displacement when t = 0 

points about which the equation is expanded in a power series 
energy loss or change 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 


AW 

n 

AW 

v 

a 

a 

n 

Oil, 0*2 

<5 

€ 

a 

T 

<P 

CO 


Definition 

energy loss for a nonlinear damping function 

energy loss for viscous damping 

damping coefficient, Case V 
th 

the n complex exponent of the exponential series 

basic exponents of the approximate solution 

energy content, AW/2W 

nonlinear spring coefficient 

exponent of e, exponential decay factor 

nondimensional time 

phase angle of displacement response 

frequency of the displacement response 


ix 



ANALYSIS OF STRUCTURAL RESPONSE WITH 
DIFFERENT FORMS OF DAMPING 

SUMMARY 


Five different analytical representations of the damping force are in- 
vestigated to provide a basis for deciding what type may exist in a given structural 
system. These representations of the damping force are incorporated along with 
a nonlinear spring force in a single-degree-of-freedom representation of the 
structural response. 

Energy loss per cycle calculations provide one basis for comparison. 

The dependence of energy losses on amplitude and frequency for sustained and 
transient response is shown. For sustained oscillations, a ratio of the energy 
loss for each damping force to that for viscous damping is plotted as a function 
of amplitude and frequency. Spring plus damping force as a function of displace- 
ment yields hysteresis loops of different shapes for different damping forces. 
Also, by selecting reasonable numerical values, transient displacement versus 
time curves and hysteresis loops are obtained by numerical integration on the 
digital computer. 

Approximate analytical solutions for each of the five cases are also shown, 
and equations for the displacement versus time response envelope are obtained. 
Considerable differences in the response envelopes are observed using the same 
numerical coefficients in each case. 

Energy loss calculations along with the displacement versus time response 
envelope and hysteresis loops provide bases for deciding what type of damping 
may exist in a structural system when observing experimental results. 


INTRODUCTION 


In the present analysis and design of missiles, the damping is generally 
considered to be viscous damping. When considering the complex system as a 
series of lumped masses, the damping force on each generalized mass is con- 
sidered proportional to velocity. When exciting the first mode of vibration, the 
modal damping is also taken as a function of velocity. As missile designs 
become more complicated, the multiple stage interfaces can produce slip 
damping, and large deflections can cause nonlinear stiffness and damping. 



The investigation of different analytical representations of damping has 
lagged far behind efforts to obtain a better representation of stiffness. To im- 
prove this situation, five different forms of damping are investigated herein. 
These five functions are each separately included in a single-degree-of-freedom 
equation with a nonlinear cubic spring force. The step input as well as the free 
oscillation condition is studied. 

Three different criteria are used to study the five different damping 
cases. Energy loss per cycle is, shown to change from case to case with ampli- 
tude and frequency. Steady state as well as transient energy losses are 
calculated. Secondly, hysteresis loops can be used to compare each case. 
Finally, numerical results are obtained by numerical integration on the digital 
computer. The envelope of each displacement versus time curve is calculated 
by first obtaining an approximate solution for the nonlinear differential equation. 

To obtain analytical representations of slip damping, it is necessary to 
match experimental hysteresis loops with those produced analytically. Thus, 
hysteresis loops for each analytical representation of the damping force shown 
herein can be correlated with current experimental efforts. 


ANALYTICAL REPRESENTATION OF DAMPING 


Consider the differential equation of motion, 


M + g , x^ + h(x) = f(t) , (1) 

where M is the generalized mass, x is the generalized mass displacement, 
g is the damping force function, h is the spring force function, and f is the 
input force. Only one type of spring force function, given by 

h(x) =K(x+ ex 3 ) , (2) 


is considered. This function is plotted in Figure i and is a typical representa- 
tion of a nonlinear spring [ 1] . In this instance, e is negative. With e = 0, 
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FIGURE 1. NONLINEAR SPRING FORCE 

the spring force is linear with a spring stiffness K. Point A on Figure 1 may 
be thought of as the point of instability in the spring force and may be used to 
select the value of e. For this study, the value € — - 1/3 is used. 

The free response and the response to a step input are considered. Thus, 


f(t) = 0 

(free) 

f(t) = p 

(step) . 


(3) 


The purpose of this study is to investigate different forms of the damping 
force function g(dx/dt, x). A study recently completed by this author 1 2] , 
discusses different sources of damping and shows different analytical repre- 
sentations of damping. The five different damping functions considered in this 
report are listed in Table I. Also shown is the possible source of each damping 
function. 
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TABLE I. ANALYTICAL REPRESENTATION OF DAMPING 


Case 


Source 

I 

Cl dx 
dt 

Viscous or 
environmental 

n 

C.lxl (sgn^) 

Material or 
slip 

HI 

C * (df )' ( Sg “ f ) 

Environmental 

IV 

/-.i , dx 

c s |xl dT 

Material 

V 


All 


To study the systems which can be described by equation ( i) , it is 
desirable to normalize the independent variable "time. " First substitute 
equation (2) into equation (1) and divide by mass M to obtain 


d 2 x ^ _i_ / dx 

diT M S \ dt 


‘) + M <X + €XS> 


P_ 

M 


for the step input. To normalize let 
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where n/k/M is the linear circular frequency. Obtaining derivatives with 
respect to t and dividing by K/M gives 


** 1 g . . o P 

X + — g(x, x) + x + ex 3 = — 


( 4 ) 


where 


x 


dx 

dr 


and x 


d 2 x 

* 


The damping force terms shown in Table I are expressed in terms of 
the normalized time as follows 


r 


^KM 


K 


I x | ( sgn x) 


^■g(x,x)= f ^(sgnx) 




n 


m 


( 5 ) 


4 KM 


I x I x 


IV 


11 ( i + G!X 2 )X , 


n/ KM 


V 
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The units of x are consistent with the other units. With mass in slugs and K 
in pounds per foot, the C's are in terms of pounds, feet, and seconds, and x 
is the displacement in feet. Any other set of units is equally applicable. 

Three different but related approaches for studying similarities of, 
and differences between, these damping terms are included in the following 
sections. 


ENERGY LOSSES AND HYSTERESIS CURVES 


A useful technique in studying different analytical representations of 
damping is to make a comparison of energy losses per cycle, AW/cycle. The 
value, AW/cycle, is obtained by integrating the spring plus damper force times 
the displacement over a vibration cycle. The spring force integrated over a 
cycle is zero; thus, 


AW 

cycle 


- J 


cycle L 


(f-) 


dx 


This equation is changed to an integral over time, or 


AW 

cycle 


= f 8 

period L 


(t- x ) 


dx 

~ dt 
dt 


( 6 ) 


Some authors prefer to use the energy content of the vibrating system, which 
is the energy loss per cycle divided by twice the vibration energy, 2W; that is, 


_ AW/ cycle 
6 ~ 2W 


AW/cycle 
dx\ 2 

max 


M 


( 


(7) 
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First consider the input force f(t) to be a type which maintains a 
constant response amplitude. The response may be approximated by 


x = A sin (cot - <f>) 


( 8 ) 


where A is the amplitude, co is the frequency, and <f> is the phase angle. The 
amplitude and frequency may be determined by using the Imethod of Klotter [3]. 
The velocity becomes 


77 = Aw cos(cot - <b) 
at 


(9) 


Substituting any of the damping functions of Table I into equation ( 6) , 
with displacement and velocity given by equations (8) and (9) , the energy loss 
per cycle is obtained by integrating between the limits of 0/co and 2ir + 0/co . 
The energy content is similarly obtained from equation (7). The results for 
each damping function are shown in Table II. 


TABLE II. ENERGY LOSSES - SINUSOIDAL RESPONSE 


Case 

AW/ cycle 

6 = AW/2W 

I 

C !7rA 2 co 

Mw 

n 

2CqA 2 

2C*. 

Mw 2 

m 

f <vA>* 

8C,A 

3M 

IV 

~ C 3 A 3 o; 
3 d 

4C 3 A 

Mco 

V 

C j7rA 2 W + ^A 4 coC 

Ciir CiiraA 2 

Mco 4Mc o 
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As shown by this table, various damping functions in effect make the energy loss 
dependent on different powers of amplitude and frequency. 

An interesting comparison of energy losses is obtained by defining an 
energy loss ratio 


R = 


AW 
n 

AW 


v 


( 10 ) 


The numerator, AW q , is the energy loss per cycle for any nonlinear damping 

function and the denominator, AW^, is the energy loss per cycle for viscous 

damping, Case I. This ratio would be the same using energy content, 6 . These 
ratios are shown in Table ID. 


TABLE IE. ENERGY LOSS RATIOS 


Case 

R = AW /AW 
n v 

R r = AW n /AW yr 

I 

1 . 0 (reference) 

measured Ci Ci 

reference Ci C . 

lr 

n 

2C 0 

2C 0 


Cprw 

C, irco 
lr 

m 

8 C 9 AC 0 

8 C 9 Aco 

3Cj7T 

3C . 7T 
lr 

IV 

4C 3 A 

4C a A 

C^ 

C lr' 

V 

1 + a f 

4 

1 +£ t ! 

4 


The commonly employed method for experimental determination of 
missile damping presupposes viscous damping, Case I. The results show different 
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values of C* as amplitude and/or frequency vary. This suggests that possibly 
some other form of damping would be more appropriate. Energy loss ratios 
provide a possible means of selecting a more appropriate damping term. 
Defining one viscous coefficient as the reference, 


Cj( reference) = , 


a reference energy loss ratio is obtained as 


AW n 

R r = AW 

vr 


where 


AW 
v 

AW 


vr 



(ii) 


These reference ratios are also shown in Table in. 

The reference energy loss ratios that are a function of frequency appear 
in Figure 2; those that are a function of amplitude are shown in Figure 3. 
Selecting one experimentally determined value of the viscous damping coefficient 
as the reference value, a plot of C^/C^ is determined as a function of both 

amplitude and frequency. Comparing this plot with those of Figures 2 and 3 
may indicate a better analytical representation of the damping coefficient. The 
coefficient would then be constant. This suggested method is to be used with 
existing experimental data and should be investigated in more detail. 

Relatively little information is available on energy losses in transient 
systems. For these systems the response can be approximated by 

x = A 0 + Ae <rt sin(ojt - $) . (12) 
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FIGURE 3. ENERGY LOSS RATIO VERSUS AMPLITUDE 
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The constant A 0 is introduced for those systems which do not vibrate about the 
null, x = 0, position. The exponential exponent, <r, must be negative for 
decaying vibrations and the final steady state (t — °°) amplitude is A 0 . A 
method for determining A 0 , A, cr, and o> in terms of the system parameters 
and initial conditions is shown in the Appendix and in the following section. The 
velocity is obtained from the derivative of equation (12) as 


~ = Acre* 7 * sin (cot - 0) + Acoe a *cos ( cot - 0) 


(13) 


The study of transient energy losses includes the following refinements 
not previously considered: 

1. The effect of a decay in the response 

2. The effect of a constant nonzero steady state amplitude 

The input forcing functions used are those given in equation ( 3) . With the free 
vibration A 0 = 0, the system is considered to be given an initial displacement and 
released. With the step input f(t) = P, the constant amplitude is considered to 
be positive. 

Substituting the various damping functions from Table I into equation (6) , 
with displacement and velocity given by equations (12) and (13) , respectively, 
the transient energy loss is obtained for each function. Care must be taken to 
integrate over the proper limits to ensure that absolute values and signs are 
correct. The results are separated into two parts. The transient energy losses 
for free response are shown in Table IV for one cycle starting at time T. The 
transient energy losses for step input are shown in Table V for one cycle 
starting at time T. 


The transient energy losses in Table IV will reduce to those in Table n 
when <j = 0. In general, the transient energy losses become increasingly im- 
portant the greater the ratio cr/co. For most applications with slow decay and 
large frequency, the assumption of a steady vibration does not yield large energy 
loss errors. However, as missiles become more flexible and slip joints at stage 
interfaces increase the rate of vibration decay, the transient energy losses should 
be considered. The transient energy losses of Tables IV and V vary with ampli- 
tude and frequency in the same way as in Table II . 
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TABLE IV. TRANSIENT ENERGY LOSSES - FREE RESPONSE 


AW/ Cycle 


_ .2 2<tT ( 7J 2tt \ w 
CiA^e le w - 11 — =— 

V ' <z) 

l 3<r or \ 

ye w 17 + e w j 


CoA 2 e 2<T ^ le 00 + e 40 


C,A 2 e 3aT * + 1 


3cr \ 2 a; 2 


■[* ♦<£)*] 




[‘ * (¥)’] 


CiA^e e w - i 


2crT / Tf27r 


Ci«A 4 ( e 40 ^ - 1 


-5 MS)*] 
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TABLE V. TRANSIENT ENERGY LOSSES - STEP RESPONSE 



I 



The transient energy losses of Table V vary considerably from those 
in Table IV when the damping term involves the displacement x. This shows 
the dependence of energy losses on the average displacement which can be 
related to average stress in many problems. 

Related to the energy loss is the force versus displacement, hysteresis 
loop, curve. In fact, the energy loss per cycle is the area enclosed by the 
hysteresis loop in the same cycle. To obtain a better understanding of the dif- 
ferent damping terms of Table I, these hysteresis loops are considered. 

The force is the spring plus damping force from equations (1) and ( 2) : 

force = g + K(x + ex 3 ) . (14) 

The sinusoidal displacement and velocity of equations ( 8) and ( 9) are substi- 
tuted into this equation, and force may be plotted as a function of displacement 
x. The result is either a single or double loop curve. Those cases that cause 
a single loop are sketched in Figure 4. Those which produce a double loop are 
sketched in Figure 5. 

The effect of slip damping is commonly obtained as a force-displacement 
hysteresis loop. A combination of curves given in Figures 4 and 5 can be used 
to approximate the slip damping condition. This in turn will allow slip damping 
to be described analytically. 

TRANSIENT RESPONSE 


Another method for comparing various analytical representations of the 
damping force is to observe displacement versus time curves and in particular 
the envelope of these curves. This requires a solution of each equation obtained 
by substituting equations ( 5) into equation ( 4) . An approximate solution is ob- 
tained in the form 


or 


x = A 0 + Ae sin(ojr - (f>) 


(15) 


which is the same form used to obtain transient energy losses (equation ( 12) ) . 
A numerical integration solution of each equation is also obtained. 
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A method for obtaining the coefficients A 0 and A, the exponent a, the 
normalized frequency w, and phase angle <p is shown in the Appendix. Two 
inputs to each system are considered, as shown in equation (3). 

For the step input, P/K, the exponent cr and the frequency oj are ob- 
tained from equations (A. 2) , (A. 6) and (A. 7) with x 10 = 0 and x 2 o = A 0 . The 
constant amplitude A 0 is obtained from equation (A. 10). The amplitude A and 
phase angle 4> are available from equation (A. 14) . Shown in Table VI are the values 
of cr and oj as well as the steady state amplitude equation for each case. 

For free oscillations, P = 0, the exponent cr and the frequency oj are 
also obtained from equations (A. 2) , (A. 6) and ( A. 7) . However, in this instance, 
the values of x 10 = 0 and x 2 o = xg/2 are found to yield the best results. The 
constant amplitude A 0 is zero. The amplitude A and phase angle <p are available 
from equation (A. 15). Table VII shows the values of cr and oj for free oscilla- 
tions resulting from an initial displacement, x 0 . 

With the approximate solution known, equation ( 15) , the envelope of the 
displacement, x, versus nondimensional time, r, curve is obtained for each case. 
This is done by defining maximum and minimum values of x and the corresponding 
times for each. Connecting these points by a continuous curve shows the envelope 
of the transient response. Maximum and minimum values are shown for the 
step input in Table VIII and for the free oscillations in Table IX. For Case n 


CO = OJ l 

A 0 = A 0 i 


( sgn x) (sgn x) = + 1 


and 


co = 0J 2 


A 0 - A^ 


(sgn x) (sgn x) = - 1 . 


For the velocity squared damping, Case III, the approximate analytical solution 
shows no decay and is not a good approximation for this case. 
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TABLE VI. SOLUTION EXPONENT, FREQUENCY, AND STEADY STATE AMPLITUDE - STEP INPUT 






to 

o 


TABLE VII. SOLUTION EXPONENT AND FREQUENCY- FREE OSCILLATION 




TABLE Vm. MAXIMUM AND MINIMUM x VERSUS T VALUES - STEP INPUT 



TABLE IX. MAXIMUM AND MINIMUM x VERSUS T VALUES - 

FREE OSCILLATION 
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Numerical results are presented in Figures 6 through 35. These results 
are presented for each of the five cases as two curves for two values of damping 
(approximately two percent and ten percent of critical) with the step input, plus 
two curves with the same two values of damping but with only an initial displace- 
ment input. Two additional curves are included for each case (with approxi- 
mately ten percent of critical damping) showing the spring plus damper force 
versus displacement, or hysteresis loops, for the same two inputs. The solid 
curves represent the fourth-order Runge-Kuta numerical integration on the 
digital computer. The dashed curves are the results of the approximate 
analytical solution shown as an envelope of the transient response. Table X is 
included as a summary of Figures 6 through 35. 

The analytical and numerical integration results for displacement versus 
nondimensional time are in good agreement. The one exception is the results 
for Case III. Horizontal deviations of maxinjum and minimum values, comparing 
analysis with numerical integration, indicate errors in the frequency to. Vertical 
deviations indicate errors in the amplitude A and the exponent <r. Considerable 
difference in the response envelope is noted when the same coefficients are used 
with different analytical representations of damping. 


CONCLUSIONS 


In this study, a single-degree-of -freedom representation of the mass 
displacement permits Ihe investigation of different analytical representations 
of damping forces. For different damping forces, the energy loss per cycle 
has different powers of amplitude and frequency. • Transient energy losses are 
considerably different than steady state losses when the vibration decay is 
large and the frequency is small. For damping forces dependent on amplitude, 
the mean response amplitude affects the energy loss per cycle. 

Different damping forces yield different hysteresis loops. It should be 
possible with a series of damping terms to approximate most experimentally 
determined hysteresis loops. Some of these hysteresis loops are dependent 
on mean vibration amplitude. 

Differences occur in the amplitude versus time curve envelope using the 
same equation except for the dependence of damping on displacement and 
velocity. In some cases, the differences are very pronounced. A three-term 
exponential series solution of the nonlinear differential equation can be used to 
obtain these curves and their envelopes. 
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It appears that a single experimental amplitude versus time envelope 
cannot be attributed to one form of damping without further investigations. In 
other words, the magnitude of the damping coefficient and the damping function 
may be simultaneously changed to yield, curves which are very similar. 

Possible procedures for determining the true analytical representation 
of damping force by one term or a series of terms are as follows: 

Energy losses could be taken as a power series in amplitude and frequency 
with coefficients determined experimentally. The correlation of existing 
viscous damping information with energy loss ratios appears feasible. 

And, finally, the actual shape of the transient amplitude versus the 
response envelope can be correlated with those generated by different 
analytical representations of damping. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsvilld, Alabama, October 17, 1966 
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FIGURE 6. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE I 

(x + 0. 04x + x - — = 0. 2) 



T =nTk/M t 


FIGURE 7. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE I 

x 3 

(x + 0. 2x + x - — = 0. 2) 



FIGURE 8. DISPLACEMENT VERSUS NONDIMENSIONA L TIME - CASE I 

x® 

(x + 0. 04x + x - — = 0, x( r = 0) = 0. 2) 

O 



FIGURE 9. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE I 

x 3 

(x + 0. 2x + x - — = 0, x( t = 0) =0. 2) 



FIGURE 10. FORCE VERSUS DISPLACEMENT - CASE I 

x 3 

(x +0.2x + x- — = 0.2) 

O 





T =Vk/M t 

FIGURE 12. DISPLACEMENT VERSUS NONDIME NSIONAL TIME - CASE II 

. x 3 

(x + 0. 04 1 x| ( sgn x) +x- — =0.2) 

O 



T = K/M t 

FIGURE 13. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE II 

x 3 

(x + 0. 2 1 x | (sgn x) + x - — = 0. 2) 

o 


I 




FIGURE 14. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE II 

x® 

(x + 0. 04 1 x | (sgn x) + x - — = 0, x(r = 0) = 0. 2) 



FIGURE 15. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE II 

x 3 

(x + 0. 2 1 x | (sgn x) + x - — = 0, x(t = 0) = 0. 2) 
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FIGURE 16. FORCE VERSUS DISPLACEMENT - CASE II 

• . x 3 

(x + 0. 2 J x | (sgn x) +x- — =0.2) 

O 



X 


FIGURE 17. FORCE VERSUS DISPLACEMENT - CASE II 

x 3 

(x + 0. 2 1 x | (sgn x) + x - — = 0, x(t = 0) = 0. 2) 





FIGURE 20. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE III 

o X 3 

(x + 0. 04k 2 + X - — = 0, x(t = 0) = 0. 2) 

o 



FIGURE 21. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE III 

x 3 

(x + 0. 2x 2 + x - -5- = 0, x(t = 0) =0. 2) 


FIGURE 22. FORCE VERSUS DISPLACEMENT - CASE III 

x 3 

(x + 0. 2x 2 + x - T- = 0. 2) 

U 



X 


FIGURE 23. FORCE VERSUS DISPLACEMENT - CASE III 

x 3 

(x + 0. 2x 2 + x - — = 0, x(t = 0) = 0. 2) 



FIGURE 24. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE IV 

x 3 

(x+0.04|x|x + x- — =0. 2) 

O 



FIGURE 25. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE IV 

. . x 3 

(x +0. 2|x|x + x - — =0. 2) 











I 



FIGURE 32. DISPLACEMENT VERSUS NONDIME NSIONAL TIME - CASE V 

x 3 

(x + 0. 039 2( 1 + x 2 )x + x - — = 0, x(t = 0) = 0. 2) 

O 



FIGURE 33. DISPLACEMENT VERSUS NONDIMENSIONAL TIME - CASE V 

x 3 

(x +0. i0(l+X 2 )x + X“~=:0, x(t=0) =0. 2) 

o 


If 
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FIGURE 34. FORCE VERSUS DISPLACEMENT - CASE IV 

„ X 3 

(x + 0. i( i + x 2 )x + x - — = 0. 2) 

u 



X 


FIGURE 35. FORCE VERSUS DISPLACEMENT 

X 3 

(x +0. i(i+x 2 )x + x- — =0, x(t=0) = 0.2)) 






APPENDIX 

APPROXIMATE SOLUTION FOR TRANSIENTS 


The approximate solution of the nonlinear differential equation of 
equation (4) is given as an exponential series 


m a t 
x = A 0 + J A n e n 
n=i 


(A. 1) 


This series when substituted into equation ( 4) yields an equation with some 
terms which are single series and others which are multiple series. Equating 
coefficients of single series terms yields a set of basic exponents cx i and Q! 2 . 
Other exponents are linear combinations of the basic exponents. The constant 
term A 0 is obtained by letting t — 00 and equating the remaining terms. 
Remaining coefficients are obtained bv evaluating m initial conditions from 
equation (4) and settingthese equal to the derivatives of equation (A.l) evaluated 
at r = 0. This method is discussed in much more detail by the author L 4 J . 
Only sufficient information is included to permit the approximate solution of 
the equations being studied. 

Essentially the method outlined is equivalent to expressing the original 
equation, equation (4) , as two first-order equations 


» 1 , . o 

— g(*i, X 2 ) - x 2 - ex§ + — 


= f 1 (x 1 , x 2 ) 


x 2 =x t 




(A. 2) 


= f 2 (x lf x 2 ) J 
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where x 2 = x and both x t and x 2 are expressible as a series. Equations (A. 2) 
can be expanded in a Taylor's series about some value x t = x 10 , x 2 = x 20 . Thus, 


x. = f.(x 10 , x 20 ) 



(xj - x 10 ) 
x 10» x 20 


+ 



(x 2 - x 20 ) + higher order terms 
x 10> x 20 

i = 1,2 . 


(A. 3) 


Since x 1 and x 2 are exponential series 


x l = a'xj = x 2 
x 2 = a"x 2 = Xj 


which are the variational equations or the linear approximations. Substituting 
the second into the first and realizing x 2 = (a') 2 x 2 yields a 1 = a". Therefore, 
in matrix notation 


x i 


a' 

0 


X 1 

— x 2 — 


_0 

a 


— x 2 - 


Expressing equation (A. 3) in matrix notation and substituting equation 
(A. 4) into this equation gives 
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0 


fl(x 10 , x 20 ) 

0 


_f 2 (xi 0 , x 20 )_ 


~dfr , 

9fj_ 



" ■ 

a ~ a 

dXf 

8x 2 



x i 

8f, 

8xj 

8f, _ 
8 x 2 

a 1 


_ x 2 . 


- 

X 10> X 20 



8fi 

- 



3xj 

8x 2 



■xio“ 

8f, 

8f? 




8x t 

8x 2 



_ x 2 0 _ 




X 10.X 20 



+ higher order terms 


(A. 5) 


Terms not involving and x 2 should be identically satisfied by the proper 
choice of x 10 and x 20 . The final steady values x 10 = 0 and x 20 = A 0 form one 
such set. Therefore, the determinate of the coefficient matrices must be zero 
to satisfy equation (A. 5). Thus, 


8fi 

8fi 

— 1- _ a t 

8xj 

8x 2 

8f, 

8f, 

8xj 

8x 2 


= 0 


o' 


x 20 
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which yields 


( o '*) 2 - O ' 1 




9f, 9f, 
9xj 9 x 2 


9fi 9f, 
9x 2 9xj 


= 0 , 


(A. 6 ) 


where partial derivatives are evaluated at x 10 , x 20 . Solving this equation yields 
two roots of a ' which are the basic exponents 


a i = cr+ ju> 

and 

a 2 = v- jco . 


(A. 7) 


The roots of d are complex when underdamped. Higher order exponents 
o? n (n > 2 ) are found to be linear combinations of a* and a 2 , as the higher 

order terms of equation (A. 5) , involve exponents 04 + oi 2 and so on. 

For the purpose of this investigation, a three-term approximation 
(m = 2 in equation (A. 1)) is used so 04 and a 2 are all that is required. Thus 


x = A 0 + Aje ** 17 + A 2 e a2T 


(A. 8 ) 


where A* and A 2 are also complex conjugate 


and 


Ai = aj + jbj 


A 2 - a t - jbj 


(A. 9) 


to ensure a real solution. 
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The coefficient A 0 is found from the steady state conditions kj = x 2 = 0 
in equation (A. 2). Thus, 


^ g(0, A 0 ) + A 0 + eA| = ^ 


(A. 10) 


becomes the steady state amplitude equation which must be solved for A 0 . The 
coefficients a t and bj are obtained using initial conditions 


x(r = 0) = x„ 

and 

x( r = 0) = x„ . 


| (A. 11) 


Therefore, substituting equations (A. 7) and (A. 9) into equation (A. 8) and 
evaluating x and x at t = 0 yields 


x 0 = A 0 + 2a t 

and 

x 0 = 2 era* - 2(obj 


That is, 



(x n - A n )o--x n 
1 ' 2cu 


and 


X = A n + 


( 


xn - An 
2 


(xn - A n ) 
2(o 


°~- Xn ^ 


(<T+jO))T 


+ 



( Xq — Aq) o’ — 
2(o 


Xfl_^ 


(o--j(o)r 

6 
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This three-term approximation may be expressed as a sinusoidal function 


err 


x = A 0 + Ae sin(ojr - (p) 


(A. 12) 


where 


A = - — [(X„ - Aq) 2 (<7 2 + CO 2 ) - 2(x 0 - A 0 )x 0 <7 + Xg ] 1 / 2 
CO 


<£ = tan 


-1 (x n - A n ) co 
(x 0 - A 0 )cr - x 0 


(A. 13) 


For the step input with zero initial conditions 


x 


= A 0 + ^ (a 2 + co 2 ) ly ^ 2 e CTT sin(coT - <p) 
co 


(A. 14) 


<p = tan 


-1 HL 


and for the free vibration with no initial velocity ( P = x 0 = 0) 


co 


(a 2 + co 2 ) 1//2 e <TT sin(cor - 0) 


(A. 15) 


4> = tan 


-1 0L 


where <r and co are the real and imaginary parts of the roots of cd from 
equation (A. 6). Phasor representations of equations (A. 14) and (A. 15) are 
shown in Figure A. 1. 
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